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ABSTRACT 


7  /-a.-;../'/ 


In  this  paper  we  consider  two-sided  parabolic  inequalities  of  the  form 
(li)  lb  1  u  <_  i//  ,  in  Q  ; 

(lii)  [  -  +  A  (t)  u  +  H(x,t,u,Du)]e  >0,  in  Q  , 

ot 


for  all  e  in  the  convex  support  cone  of  the  solution  given  by 
K  (u)  =  {A(v  -  u)  :  <_v  £  ,  X  >  0}  ; 


(li ii) 


0 


u ( •  ,T)  =  u 


where 

Q  =  Q  x  (0,  T)  ,  l  =  x  (0,T)  . 

Such  inequalities  arise  in  the  characterization  of  saddle-point  payoffs  u  in 
two  person  differential  games  with  stopping  times  as  strategies.  In  this 
case,  H  is  the  Hamiltonian  in  the  formulation.  A  numerical  scheme  for 
approximating  u  is  obtained  by  the  continuous  time,  piecewise  linear, 
Galerkin  approximation  of  a  so-called  penalized  equation.  A  rate  of  conver¬ 
gence  to  u  of  order  Oth'*"^)  is  demonstrated  in  the  L^(0,T  ;  H'*’  (.';))  norm, 
where  h  is  the  maximum  diameter  of  a  given  triangulation . 
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SIGNIFICANCE  AND  EXPLANATION 


\  ! 
i 

""^This  paper  studies  the  two-sided  differential  inequality  which  charac¬ 
terizes  the  optimal  payoff  in  a  two-person  differential  game,  when  the 
strategies  available  to  the  players  are  stopping  times.  By  converting  the 
inequality  to  an  equation  via  a  so-called  penalization  method,  we  are  able 
to  apply  a  standard  numerical  method  for  approximating  solutions  of  nonlinear 
equations.  We  obtain  new  rates  of  convergence  for  the  procedure 
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CONVERGENT  APPROXIMATIONS  IN  PARABOLIC  VARIATIONAL 
INEQUALITIES  II:  HAMILTON- JACOBI  INEQUALITIES 


Joseph  W.  Jerome 


§1 .  Introduction. 

In  this  note  we  continue  the  investigation  of  the  discretization  of  parabolic 
variational  inequalities  begun  in  [3].  The  problem  is  decidedly  more  general  here 
and  our  approach  correspondingly  different.  As  distinct  from  the  one-sided 
inequalities  defined  by  linear  differential  operators  such  as  the  heat  operator, 
studied  in  (3],  we  consider  here  two-sided  inequalities  defined  by  nonlinear 
operators.  Thus  we  seek  a  function  u  satisfying,  on  a  space-time  domain 
Q  =  x  (0  ,T)  : 

(i)  ij/^  <  u  £  \li2  ,  in  Q  ; 


H  ( x  ,  t ,  u  ,  Du )  =  0  , 

if  ^  <  u  <  'v2  ; 

H (x, t , u ,Du)  >  0  , 

if  u  =  i(i  ; 

H  (x,  t  ,u,Du)  0  , 

if  u  =  <p2  ■ 

Adjoined  to  (l.li,  ii,  iii,  iv)  is  a  standard  homogeneous  Neumann  boundary  condi¬ 
tion  and  a  terminal  condition  at  t  =  T  : 

(l.lv)  =  0  on  J  =  x  (0 , T)  ;  u(>,T)  =  u  . 

3v  L 

The  precise  hypotheses  are  stated  in  section  two,  where  we  define  a  penalized 

1/4 

problem  with  solution  u^  and  demonstrate  a  rate  of  convergence  of  order  0(c  ) 

2  1 

in  the  L  (0,T  ;  H  (SI))  norm  (cf.  Theorem  2.2).  In  section  three  we  define  the 
Faedo-Galerkin  approximation  of  u£  and  demonstrate  convergence  to  u  with 

the  rate  of  0(h1/^^)  in  the  L^(0,T  ;  H^(SD)  norm,  if  piecewise  linear  elements 

3/4 

of  maximal  diameter  h  are  employed  with  h  =  e  (cf.  Theorem  3.3). 
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(ii)  -  +  A (t) u  + 

ot 


(1.1) 


(iii) 


f  +  A(t,u  + 


(iv)  -  +  A ( t )u  + 


The  formulation  (1.1)  arises  in  stochastic  control  and  two-player  differential 
game  theory,  where  it  is  realized  as  a  (double)  minimax  value  of  a  stochastic 
functional  (cf.  Bensoussan  and  Lions  12)).  In  contradistinction  to  the  Hamilton- 
Jacobi  equation,  which  arises  as  a  (single)  minimax  value  associated  with  the 
Legendre  transformation,  the  inequality  (1.1)  involves  the  introduction  of  an 
optimal  stoppinq  time.  For  the  sake  of  brevity,  we  omit  the  details  and  instead 
refer  the  reader  to  (21.  Nonetheless,  we  should  observe  the  nonstochastic  applica¬ 
tions  of  (1.1),  such  as  certain  Ftofan  problems  (cf.  Bensoussan  and  Friedman  [lj). 
The  methods  of  this  paper  are  obviously  applicable  to  the  simpler  model  discussed 
in  13].  We  have  chosen  a  (natural)  Neumann  boundary  condition  rather  than  a 
(forced)  Dirichlet  boundary  condition  to  simplify  the  exposition  of  / 3 ,  since 
piecewise  linear  elements  cannot  vanish  on  an  arbitrary  surface  •  .  Wo  mention 
finally  that  higher  rates  of  convergence  are  to  Ik-  expected  if  U  'b  viz., 

Op  1^2)  convergence  of  the  solutions  of  the  serialized  o  stations  and  o(h  ) 
convergence  of  tin-  finite*  element  »;?roXHnut  ions,  with  h  ~  •  • 


§2 .  The  Penalization. 

N 

Let  ft  be  a  smoothly  bounded  domain  in  3R  and  let  A(t)  be  a  uniformly 
coercive,  elliptic  operator  in  divergence  form: 

(2.1)  A  ( t)  =  -  7  D.  a.  .  ( •  1 1)  D .  +  V  b,  (.,t)D,  +  c(-,t) 

.  .in  i  L  k  k 

i ,  3  J  k 

where  D  =  ~  ,  a.  .  £  C1(Q)  and  b  ,  c  e  C1([0,T]  ;  C(ft))  and 

i  ox.  k 


2 

(2.2)  a  (u.u)  >  c  ||  u || 

H  (ft) 


for  all  u  e  H1 (ft) .  Here 


cQ  >  0  , 


(2.3) 


a  (u,u) 


i.D 


(D.U, 


ID 


D.u) 

3  L2(ft) 


+  E  (bk  V’U)  2 

k  (ft) 


(cu,u)  , 

L  (ft) 


-  N-t-1 

Let  H  :  ft  x  [0,T]  x  K  -+  3<  bt  a  continuous  function  satisfying,  generically, 
(i)  |h(x,s,v,Dv)  -  H(x,t,w,Dw)j  £C(|s  -  t|  +  |v  -  w|  +  |dv  -  Dw|j  , 

(2.4) 

(ii)  |  H  (x,  t,  v,Dv)  |  <_h(x,t)  +  |  v  |  +  |  Dv  |  , 


where  D  =  (D, ,  ...  ,  D„)  and  h(’,.)  is  a  bounded  measurable  function;  I* I 
IN 

denotes  an  appropriate  Euclidean  norm.  Let  ^  and  ^  be  given  satisfying 


(2.5) 


(i)  i/l  £  H2  (Q)  ,  i  =  1,  2, 

(ii)  l/'1  £  '  in  2  • 


We  shall  define  a  class  of  penalized  problems  depending  on  a  parameter  E  >  0. 
These  provide  both  a  tool  for  proving  the  existence  of  a  solution  of  (1.1)  in  the 
class 

(2.6)  Xx  =  L2(0,T  ;  H2 (ft) )  n  H1 [0,T  ;  L2(ft)] 


as  well  as  defining  the  base  equation  for  the  Faedo-Galerkin  approximation  of  the 
next  section.  This  penalized  equation  is  exactly  the  one  introduced  in  [21. 
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We  define  now  the  penalized  problem,  for  e  >  0, 


l  +  i  _ 

( i) - +  A (t) u  +  — (u  -  lit.) - (u  -  )  +  H(x,t,u  ,Du  )  =  0  ,  in 

9t  £  £  £  2  £  £  1  E  £ 

3u 

(2.7)  (ii)  — ^  =  0  ,  on  l  , 

(iii)  u  ( *,T)  =  u  , 
e 

The  functions  (*)+  and  (•)”  are  the  positive  and  negative  parts  of  the  identity, 
defined  so  that 

+  —  1  4-  — 

t  =  t  -  t  .tern  ;  t  ,  t  >_  0 

Remark  2.1.  It  is  known  that  (2.6)  possesses  a  solution  u  satisfying 
-  £ 

(2.8)  u  €  L2(0,T  ;  H2(M)  n  H2[0,T  ;  L2  (fi)  ]  =  X 

£  2 

if  say 

2 

(2 . 9i )  u  £  H  (Q)  . 

If,  in  addition, 

(2.9ii)  i|i  (*,T)  <_  u  <_  i|i2  ( - , T)  , 

then  one  can  use  the  equations  (2.7)  to  prove  the  existence  of  solutions  u  of 
(1.1)  in  the  regularity  class  of  (2.6),  under  the  hypotheses  described  earlier 

in  this  section  by  the  arguments  of  the  proof  of  Th^or&me  1.1,  p.  449,  of  [2] 

(cf.  [2,  pp.  449-455]).  Actually,  one  can  prove  the  existence  of  u  £  X^  satisfy¬ 
ing  (l.li,v)  and 

T 

(2.10)  /  (A(t)u  +  H ( • , t , u , Du )  -  It  .  v  -  u)  dt  >  0 

0  I .  <SJ> 

2 

for  all  v  f  L  (p),  <_  v  <_  i ,  from  which  (1 .  Iii, iii ,  iv)  follow  if  and  g^ 

are  coincident  only  on  a  set  of  measure  zero: 

(2.11)  meas  { i|i ^  =  i)^}  =  0  (in  Q) 

Standard  methods  also  give  the  characterization  of  (1)  of  the  abstract. 


Remark  2.2.  As  remarked  in  [2] ,  there  is  no  loss  of  generality  in  assuming  that 

12  2 

H(*,  •«•,♦)  defines  a  monotone  mapping  from  H  (£2)  c  L  (£2)  into  L  (£2)  : 


(2.12) 


(H (  •,  ■ , v ,Dv)  -  H(  *,  •  ,w, Dw)  ,  v  -  w)  ^  0 

L  ( £2) 


for  all  v,  w  e  H  (£2).  We  shall  assume  (2. 1-2. 5)  and  (2.8,  2.9,  2.11,  2.12)  for 
the  sequel.  The  reduction  to  (2.12)  is  achieved  by  an  integrating  factor  and  change 
of  variable. 

Proposition  2.1.  The  following  'a  priori'  estimates  hold  for  the  solutions  of  (2.7): 
3u 


(i)  || 


L  (Q) 


+  II  «  II  2  i  £  C  ' 

L  (0  ,T  ;  H  (£2)  ) 


(2.13) 


<ii>  ||  «.  II  ,  2  i  C 

L  (0,T  i  H  (£2)) 


(iii)  ||  (u£  -  *2)  |Il2(q)  +  II  (»c  "  *!> 


L2(Q> 


<  c 


(iv)  7f  11  £  <ue  -  V  +  II  2  +  ji  II  ,ut  ‘  V' 

L  / 


L2(2) 


<  c 


Proof:  The  estimates  (2.1 3i ,  iii,  iv)  are  derived  as  in  (2,  pp.  449-455] .  Estimate 
(2.13ii)  follows  from  these  and  (2.7). 

Remark  2.3.  These  estimates  are  sufficient  to  prove  the  existence  of  a  solution 

32u 

u  of  (1.1)  in  X  .  Note  that  the  existence  of  - &  was  used  to  obtain  the 

.  3t 

3u 

’a  priori'  estimate  of  — ~  . 

dt 

Theorem  2.2.  The  solutions  of  the  penalized  problems  (2.7)  converge  to  the  solution 
of  (1.1)  with  order  0(e^'S  in  the  norm  of  L2(0,T  ;  H^(£2))  : 


(2.14) 


T  2 

/  II  u  -  UF  II  i  dt  <  C  c 
0  H  (£2) 


1/2 


Proof:  We  note  that  the  function  r^_,  defined  by 
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1 


(2.15i) 
satisfies 
(2 . 15ii) 


r£  =  uc  -  <ue  -  *2)  +  (tt£  - 


h  -rc-  *2 


Thus,  from  (2.10)  we  have,  using  integration  by  parts. 


(2.16) 


.  3u 

+  H(-,t,u,Du)  ,  r  -  u)  +  a  (u  ,  r  -  u)  >  0 

L2(fi)  tee 


for  almost  all  t,  0  <  t  <  T.  Multiplying  (2.7)  by  r^  -  u  and  subtracting  (2.16) 
from  the  resultant  integrated  equation  gives 
Id  2 

~  2"  dt  !1  u  “  ue  II  2  +  afc(u  -  u  ,  u  -  u  ) 

L  (ft)  Z  E  E 

+  ( H ( • , t , u  , Du  >  -  H ( • , t , u , Du )  ,  u  -  u) 

E  L2((2) 

<  (A(t)u  -  J~  ,  (u  -  ^  )  -  (u  -  it  )+)  +  (H ( • , t,u  , Du  )  -  H ( • , t , u , Du )  , 

.  l/(Q)  C  £ 

(u  -  it  )  -  (u  -  it,)+) 

E  2  L2(fJ) 


(2.17) 


where  we  have  used  (cf.  [2,  p.  2091) 


-  (u  -  it.)  •  (r  -  u)  >  0  ,  (u  -  1 1,)+  •  (r  -  u)  >  0  . 

c  1  e  g  2  e  — 

Integration  of  (2.17)  over  (0,T),  together  with  (l.lv),  (2.4i),  (2.7iii),  (2.12) 

and  the  elementary  inequality, 

(2.18) 

yields 


!<'.<»>  I  i  J  <h||  f  II  2  ♦  n'1  tl  g II  2>  , 


/  at(u  -  u^  ,  u  -  u£)dt  <C|  ||(u£  -  ltd'll  2  +  |'(u£  -  ita)  1 1  2 


l  (g) 


L  (0) 


'  a. 


'“t  •  V  2 

L  (0) 


and  the  proof  of  (2.14)  is  concluded  by  use  of  (2.1 3i ii ) . 
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§3 .  Continuous  Time  Finite  Element  Approximations. 

For  h  >  0,  let  T  be  a  triangulation  of  the  given  domain  fi.  Thus, 


(3.1) 


D  =  U  T 

tel 


where  t  is  a  typical  (closed)  element  in  the  simplicial  decomposition  T  ;  in 

h 

particular,  we  permit  nonsimplicial  elements  near  the  boundary.  Let  M,  denote 

h 

the  linear  space  of  continuous  piecewise  linear  trial  functions  determined  by  : 


(3.2) 


=  (X  <r  C(fi)  :  X| 


is  linear  Vt  €  T 


Let  E  be  the  Ritz-Galerkin  H^(S!)  projection  defined  by 


( 3 . 3  i ) 


(E  g? ,  X)  =  (*,X> 

H  (O)  H  (Q) 


for  all  X  c  M 


.1  . 


for  each  fixed  g  c  H  (T;) ;  here  we  use 
(3. 3ii) 


(g,ii)  ,  =  (Vg,Vi(i)  +  j—  /  <f  J  <Ij  , 

h  (si)  l  (q)  a  si 


which  defines  a  norm  equivalent  to  the  standard  H  (Si)  norm  in  the  usual  way 
(cf.  (4)).  Let  I  denote  the  interpolation  operator.  We  shall  assume: 


(3.4) 


*1  -  Zh  *1  -  \  *2  -  *2 


Roughly  speaking,  (3.4)  asserts  that  is  smooth  and  convex  and  is  smooth 

and  concave;  indeed,  these  assumptions  guarantee  (3.4).  We  make  the  standard 
finite  element  assumptions  (cf.  [5]) 


(3.5) 


H3  (SI) 


ch2‘j  |M 


h2(Q) 


j  =  0,  1,  £  H  (fi)  , 


for  F.  =  K  and  for  F,  =  I, 
h  h  h  h 


We  are  now  in  a  position  to  define  the  finite  element  approximation  via  a 


standard  Faedo-Galerkin  method  based  upon  (2.7). 


Definition  3.1.  The  finite  element  approximation  u  :  [0,T]  ■»  M,  is  the  uni'-iue 

-  e  h 

element  in 

X  =  H1[0,T  ;  L2[Q)]  n  L2(0,T  ;  H1(o)) 


satisfying 

(3.6i) 


<J  (u  )  ,  X>  =  0  ,  for  all  X  e  M  , 

et  h 


(3 . 6 i i ) 


u  ( • ,  T)  =  I  u 
£  h 


Here,  J  is  the  map, 
£ 


J  :  X  -*■  L2  (0 ,  T  ;  H_1  ( £7)  ) 
E  1 


defined  by  the  pointwise  relation  on  (0,T)  : 

(3.6iii)  <J  (<£)  ,  =  (-  +  G(e,^,iJ)  )  ,  i)j)  ,  +  a  =  0  , 

e  3t  12  l2<si>  1 


where 
(3 . 6iv) 


G(c  -  —  g  (ifi  ,i^  ,\p2)  +  H  ( ■ ,  •  ,\f  ,D  £)  , 


( 3 . 6v ) 


=  “fr?  -  il^)  +  (if  ~  ip2> 


Remark  3.1.  The  existence  of  a  unique  solution  of  (3.6)  is  standard  and  may  be 
achieved  via  the  theory  of  pseudomonotone  operators. 


Proposition  3.1.  The  following  'a  priori'  estimates  hold  for  the  solutions 

{uht  of  (3.5): 

C  £ ,  h 

3uh 

(i)  ll-jrll  2  i  C1  h  /E  +  c,  • 

3t  L2(Q)  1  2 

(3.7)  (ii)  II !|  2  ,  1C  , 

L  (0 , T  ;  H  (Q)  ) 

(iii)  jp  ||g(Up,<i  ,<1>  )  ||  ,  1C  , 

/c  1  2  L  (Q) 

where  C,  C^,  are  independent  of  c  and  h. 
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Proof:  Select  X  =  u  -  I,  Ji,  in  (3.6i).  Then,  we  have 
e  hi 


(3.8) 


h»;h22  +.t(u;.»5>  iu«;- v'h22 

L  (IJ)  L  C:) 

i  (-^ue  •  w  > +  vuc  -  rh  v  +  (h(u£>  '  Jh  v  2r  i 

L  (St;  l-»  •  / 


where  we  have  used  (3.4)  after  an  addition  and  subtraction,  have  written  H(u  ) 
for  the  last  term  in  (3.6iv)  and  have  noted 

<ue  '  Th  *1’  ‘  (uc  '  V+i°  ‘ 

Now  integrate  over  (t,T) .  Integrating  the  first  term  on  the  r.h.s.  of  (3.8)  by 
parts,  applying  ( 2 . 4 i i )  to  the  third  term  and  estimating  the  first,  the  second 
and  part  of  the  third  term  by  the  inequality  (2.18)  we  obtain,  for  appropriate 
choices  of  n  , 


1^1*, t)  ||  ,  +  /  (at(u^  ,  uh£)  +  ~  ||(u*  -  v1)  ||  2  } 


dt 


L  0!)  t 


L  (Q) 


(3.9) 


T  h  2 
h  M 


i  C  +  C  /  ||u  ||  2  dt  . 

1  t  E  L  (ft) 


Gronwall 1 s  inequality  applied  to  (3.9)  yields,  in  particular 
,  2  , 2 


(3.10) 


“lla  1  +  £  11  (u£  -  v  11  2  i c  • 

L  (0,T  ;  H1  (Q) )  L  (Q) 


A  parallel  argument,  with  X  =  u^  -  1^  4^ '  yields 


(3.11) 


\\  ll2  !  +  F  l|lU£  -  V+!l  2  • 

E  L  (0,T  ;  H  (ft))  L  (Q> 


Clearly,  it  remains  only  to  show  that 

„3>  h 

(3.12) 
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3u 


3^- 


Setting  X  s 
integrating  over  (0,T)  , 


r-—  in  (3.6i)  we  have,  after  adding  and  subtracting  -r-~  and 
dt  at 


3uh  2  2 

|-“  |J  +  a(u  ( •  /0)  /U  ( *  ,0)  )  +  -  I)  (uV,0)  -  (-,0))"  || 

l/(0)  E  E  E  E  1  b  (fl) 


(3.13) 


„  h 

1  -V,  h  3u. 

<  -  ((I.  u  -  <1  (•,  T) )  ,  -r-r-  )  _  +  a  (u,u)  +  (H(u  )  ,  -r^-  )  , 

~E  h  1  3t  L2(:>)  e  3t  L2(Q) 


However,  by  (2.9ii), 


a,  u  -  -  (i.  u  -  i.  e,  (%t)>  +  (i  («,t>  -  e,(*/T)) 

hi  hhl  hi  1 


-  (i  *_  ( • , t)  -  e,  ( • » t)  ) 

hi  1 


Thus,  by  (2.5i),  (3.5),  and  (2.18)  we  obtain  from  (3.13), 

u2 


(3.14) 


3uh  2 

-nr  II  ,  1  C  (  —  +  a  (u » u)  +  C)  , 

3t  L(O)  1  E 


from  which  (3.12)  follows  and  the  proof  is  concluded. 

Corollary  3.2.  There  is  a  constant  C,  independent  of  e  and  h,  such  that 
(3.15) 


||j  (uh)  ||  ^  C[h3/c  +  h /ft  +  h) 

C  E  F 


where  F  denotes  the  dual  of  H  (ft) . 


Proof:  Let  v  r  h  (ft).  Then,  by  (3.3),  we  have 


(J  ( u ^)  ,  v)  =  (J  (uh)  ,  v  -  E,  v) 
c  e  e  e  h 


3u 


=  <_  IT  +  G<*-VW  ,  V  -  Eh  V) 


1 


so  t  ha  t 


(3.16) 


J  (u  )  II  =  sup 
G  E  F  ||v  ||  ,  =1 

h  (n) 


|<  J  (u  )  ,  v> | 
e  e 


<_  C  (h  /e  +  l//e  +  C  )  sup 

”  1  2  I.  || 


I  1 
H  (fi) 


=  1 


|v-Eh  v 


Now  (cf.  [5]) 


sup 
llv  II 


llv  ~  E  v  ||  <  Ch 

=1  L  (ft) 


so  that  (3.15)  is  immediate  from  (3.16). 

Remark  3.2.  We  are  now  prepared  to  state  the  major  result  of  the  paper. 
Theorem  3.3.  There  is  a  constant  C,  independent  of  e  and  h,  such  that 


1/2 


(3.17) 


|U£  "  Uc  11  2 

L  (0,T  ;  H  (ft)) 


!  <  CIC/  ||j  <u^>  ||  dt)  ]  . 


In  particular, 
(3 . 1 8 i ) 


|u  -  u  ||  £  C(h  /e  +  h/v'e  +  /e  J 

E  L  (0,T  ;  H  (ft)) 


If  the  choice  h  =  e 
(3 . 18ii) 


3/4 


is  made  then 

I  h  ||  1/3 

lu  -  ur  II  2  i  £  C  h 

L  (0,T  ;  H  (ft)) 


Proof:  In  light  of  Theorem  2.2  and  Corollary  3.2,  it  suffices  to  prove  (3.17) 

Using  the  monotonicity  of  g(  •  and  (2.12)  we  have, 

.  h  h,  ,  ,  .  „  ,  h,  h. 

a(u  -  u  ,u  -  u  )  <(J  (u  )  -  J  u  ,  u  -  u  > 

te  e  e  c  —  c  c  ec  c  c 

=  <  J  <uh)  ,  u  -  uh  > 

E  E  EE 

<  ||J  (uMl  ||u  -  ii*  || 

F  E  £  H  (Si) 

so  that  (3.17)  follows. 


>11- 


Remark  3.3.  If  H  5  0,  then  (2.13iii)  can  be  strengthened  so  that  /T  is 
replaced  by  e.  This  leads  to  a  rate  of  convergence  in  Theorem  2.2  of  order 
(cf.  the  proof  in  the  case  of  one  obstacle  in  [2,  p.  224)).  The  'a  priori' 
estimates  of  §3  remain  unchanged  and  the  choice  e  =  h  leads  to  a  convergence 
rate  of  /h,  replacing  /h  in  (3.18ii). 
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